Anomalous Hall effect in a two-dimensional electron gas 
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The anomalous Hall effect in a magnetic two-dimensional electron gas with Rashba spin-orbit 
coupling is studied within the Kubo-Streda formalism in the presence of pointlike potential impu- 
rities. We find that all contributions to the anomalous Hall conductivity vanish to leading order in 
disorder strength when both chiral subbands are occupied. In the situation that only the majority 
subband is occupied, all terms are finite in the weak scattering limit and the total anomalous Hall 
conductivity is dominated by skew scattering. We compare our results to previous treatments and 
resolve some of the discrepancies present in the literature. 
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I. INTRODUCTION 

In 1879, Edwin Hall ran a current through a gold foil 
and discovered that a transverse voltage was induced 
when the film was exposed to a perpendicular magnetic 
fieldii The ratio of this Hall voltage to the current density 
is the Hall resistivity. For paramagnetic materials, the 
Hall resistivity is proportional to the applied magnetic 
field and Hall measurements give information about the 
concentration of free carriers and determine whether they 
are holes or electrons. Magnetic films exhibit both this 
ordinary Hall response and an extraordinary or anoma- 
lous Hall response that does not disappear at zero mag- 
netic field and is proportional to the internal magneti- 
zation: i?Haii = RqH -\- RgM, where i?Haii is the Hall 
resistance, Ro and Rs are the ordinary and anomalous 
Hall coefficients, M is the magnetization, and H is the 
applied magnetic field. The anomalous Hall effect (AHE) 
is the consequence of spin-orbit coupling and allows an 
indirect measurement of the internal magnetization. 

Despite the simplicity of the experiment, the theo- 
retical basis of the AHE is still hotly debated and a 
source of conflicting reports.-^ Different mechanisms con- 
tribute to the AHE: an intrinsic mechanism and extrinsic 
mechanisms such as skew-scattering and side-jump con- 
tributions. The intrinsic mechanism is based solely on 
the topological properties of the Bloch states originating 
from the spin-orbit-coupled electronic structure as first 
suggested by Karplus and Luttinger.- Their approach 
gives an anomalous Hall coefficient Rs proportional to 
the square of the ordinary resistivity, since the intrin- 
sic AHE itself is insensitive to impurities. The skew- 
scattering mechanism, as first proposed by SmitA^ relies 
on an asymmetric scattering of the conduction electrons 
by impurities present in the material. Not surprisingly. 



this skew scattering contribution to Rs is sensitive to the 
type and range of the scattering potential and, in contrast 
to the intrinsic mechanism, scales linearly with the diag- 
onal resistivity. The presence of impurities also leads to 
a side-step type of scattering, which contributes to a net 
current perpendicular to the initial momentum. This is 
the so-called side-jump contribution, whose semi-classical 
interpretation was pointed out by Berger.— However, it is 
not trivial to correctly account for such contributions in 
the semiclassical procedure, making a connection to the 
microscopic approach very desirable. 

The early theories of the AHE involved complex cal- 
culations with results that where not easy to interpret 
and often contradicting each other. ^ The adversity fac- 
ing these theories stems from the origin of the AHE: 
it appears due to the interband coherence and not just 
due to simple changes in the occupation of Bloch states, 
as was recognized in the early works of Luttinger and 
Kohni^ Nowadays, most treatments of the AHE either 
use the semiclassical Boltzmann transport theory or the 
diagrammatic approach based on the Kubo-Streda linear- 
response formalism. The equivalence of these two meth- 
ods for the two-dimensional Dirac-band graphene system 
has recently been shown by Sinitsyn et aZ.,12, who explic- 
itly identified various diagrams of the more systematic 
Kubo-Streda treatment with the physically more trans- 
parent terms of the semiclassical Boltzmann approach. 

It is therefore important to also obtain a similarly co- 
hesive understanding of the AHE in other systems such 
as the two-dimensional (2D) spin-polarized electron gas 
with Rashba spin-orbit interaction in the presence of 
pointlike potential impurities, where a series of previous 
studies has led to a multitude of results with discrepan- 
cies arising from the focus on different limits and/or sub- 
tle missteps in the calculations i^^'^^d^d'^ds, 16^17 jg ^j^g 
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purpose of this paper to review and analyze the previous 
attempts and to provide a detailed analysis of all contri- 
butions to the AHE in a two-dimensional electron gas. 
Since we have already demonstrated the equivalence of 
the Kubo-Streda formalism and the semiclassical Boltz- 
mann approach with respect to skew scattering in the 
two-dimensional electron gas in a previous paper, we 
will focus here exclusively on the diagrammatic formal- 
ism based on the Kubo-Streda treatment. 

The outline of the article is as follows. We start by re- 
viewing and commenting on previous studies of the AHE 
in the two-dimensional electron gas in Sec. [TTl where we 
compare them with our results and discuss the discrep- 
ancies and their possible origins. In Sec. IIIII we present 
details of our calculation within the diagrammatic Kubo- 
Streda formalism. In Sec. lIII O we provide simple analyt- 
ical limits of all terms of the anomalous Hall conductivity 
and discuss the full evaluation in Sec. IIII Dl Finally, in 
Sec. II VI we present our conclusions. 



II. COMPARISON WITH PREVIOUS 
APPROACHES 

Currently there are several publications on the AHE in 
two dimensional systems reaching different quantitative 
predictions even in the same limits . "'^"'"'^"'^i-'^^'-'^'^i-'^^'-'^^i"'^^'"'^'' In 
the present paper we present a calculation with conclu- 
sions that are in disagreement with some previous stud- 
ies. On such a background we believe that previous ar- 
ticles have to be discussed in some details. Below we 
review the history of the problem and explain why we 
think the subject has to be reconsidered. 

A first study of the AHE in two dimensional systems 
was done by Culcer et al.^ who calculated only the in- 
trinsic contribution to the Hall conductivity for a wide 
class of two-dimensional systems, including the Rashba 
two-dimensional electron gas as a special case. The in- 
trinsic contribution plays a special role in the theory of 
the AHE because it is not related to the scattering of elec- 
trons but is rather caused by the unusual trajectories of 
electrons under the action of the electric field. However, 
the disorder contributions can also be important and fur- 
ther insight was needed in the quest for a quantitatively 
rigorous theory of the dc-AHE. 

The first attempts to understand the disorder effects 
where done independently by two groups jiSiii^ each em- 
ploying different approaches. Dugaev et. ali^ used the 
version of the Kubo formula, which expresses the Hall 
conductivity in terms of the causal Green functions. The 
intrinsic contribution appears as a result of calculations 
with bare Green functions, while disorder effects renor- 
malize the quasi-particle life time and the current vertex. 
This approach is formally rigorous and is similar to the 
one we adopt in our work. However, our final results are 
quantitatively different from those found in Ref. [HI due 
to a subtlety in the calculation of the vertex at the Fermi 
surface which was later corrected in the appendix of Ref. 



[Tol . Starting with the equation for the renormalized ver- 
tex Tx = akx + bax + cay and with the assumption that 
the density of impurities is low, they find correctly that 
6 = to leading order in rn, i.e. a/6 oc n^. However, such 
a term gets multiplied by an equivalent divergent term 
within the Kubo formula leading to a non-zero contribu- 
tion to the AHE conductivity to zeroth order in n^. 

In contrast to the previous quantum mechanical ap- 
proach, Sinitsyn et alJ^ employed the semiclassical wave- 
packet approach focusing only on the understanding of 
the side-jump contribution and formulating the semi- 
classical problem in a gauge invariant form. This worki^ 
intentionally avoids a discussion of the skew-scattering 
contribution due to the asymmetry of the collision term 
kernel, which is also an important mechanism of the Hall 
current and can even be parametrically similar to all 
other contribution^" in the case of Gaussian correlations. 
Therefore, the work in Ref. [l3| is meant as an intuitive 
introduction into the physics of the anomalous velocity 
and the side-jump effect, but does not offer a rigorous 
quantitative comparison even in the considered limit of 
smooth disorder potential. 

Subsequently two papers by Liu et aZj ^'^i^^ studied the 
problem using the Keldysh technique for linear transport. 
The Keldysh technique leads to the quantum Boltzmann 
equation for the diagonal elements of the density matrix 
in momentum space when only elastic scattering events 
are considered. In the steady state limit of a weak electric 
field this equation can be written as follows: 

eE-Vp/5(p) +?:[Ho,/5(p)] = Icoi{p{p)) , (1) 

where Icoi contains all disorder dependent terms that be- 
come zero when p{p) is the density matrix in thermo- 
dynamic equilibrium and Hq is the disorder free part of 
the Hamiltonian. The "hat" means that p and Icoi are 
matrices in the band index space. The term containing 
the electric field is called the driving term. In the linear- 
response approximation it only depends on the equilib- 
rium part of the density matrix. 

To start with Eq. ([T]) is correct and is also the starting 
point of the pioneering work by Luttinger— and there- 
fore one can compare it directly with steps taken by 
Liu et ali^'^ Luttinger's approach was to split the den- 
sity matrix into equilibrium and nonequilibrium parts 
P — Pcq + Pneq whcrc p^eq IS linear in electric field. It 
is this part of the density matrix that is responsible for 
nonzero currents. For weak disorder potential V ^ Lut- 
tinger looked for p^eq as a series in powers of the disorder 
potential. He found that this series starts from the term 
of the order V"'^ 

Pneq — Pneq ' Pneq ' Pneq ^ l^j 

As pointed out by Luttinger, the leading order term pneq 
does not contribute to the Hall effect and is only responsi- 
ble for the longitudinal diffusive current. The term pneq' 
was identified with skew scattering. This term, however. 
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is parametrically very distinct and vanishes in the ap- 
proximation of purely Gaussian correlations of disorder 
Fourier components; therefore, Luttinger went to next 
order and calculated the term pneq- He found a num- 
ber of contributions, whose physical meaning he did not 
clarify. The main conclusion was that at this order both 
the diagonal and off-diagonal parts of the density matrix 
become nonzero and contribute to the Hall conductivity, 
which becomes formally independent on the strength of 
disorder V in the DC limit, although disorder has to be 
included in the intermediate calculations. 

Comparing this with the first work of Liu and Leii^ we 
find that they determined self-consistently only the off- 
diagonal part of the density matrix in band index. This 
is, however, not enough for a rigorous quantitative result 
because the diagonal part of the pnlq contribution has 
been known to be important since Luttinger's pioneering 
work. 

In their next effort Liu et al^ studied the problem of 
2D Rashba systems in small gap semiconductor materi- 
als, in which a projection to the conduction band leads to 
extrinsic type spin-dependent contributions. In this work 
they noticed that the diagonal part is important and cal- 
culated it numerically. For the driving term in Eq. ([T]) 
Liu et al. assume that peq is just a diagonal equilibrium 
Fermi distribution. This would be correct if one was us- 
ing the basis of the eigenstates of the full Hamiltonian 
with impurities. However, both Liu et al. and Luttinger 
work in the chiral basis of the disorder free Hamiltonian 
Ho- In this basis the equilibrium state density matrix is 
no longer diagonal and can also be written as a series in 
powers of the disorder potential: 

p,^ = pi^)+p(^) + ... (3) 

Luttinger has shown that in order to properly evaluate 
the non-equilibrium part p^neq one should include the sec- 
ond term peq of the expansion of the equilibrium den- 
sity matrix in Eq. ^ into the driving term of Eq. ([1]). 
This was not done in Ref. [HI and therefore we believe 
that their work is incomplete due to such omission. We 
also note that the correction of order in Eq. ^ leads 
to the Hall current contribution, which was identified in 
the semiclassical approach^^. as the anomalous distribu- 
tion correction and if omitted leads to errors of factors of 
two in the typical side-jump type contributions.^ In the 
Kubo formula approach, neglecting this correction would 
be equivalent to the unjustified omission of an important 
subset of Feynman diagrams.^" Within the calculation 
presented here all these terms are present. 

Inoue et al}^ calculated the AHE contribution using 
the same approach we use focusing on the limit of both 
subbands being occupied and, in addition to the disor- 
der that we consider, incorporating magnetic impurities 
in the model Hamiltonian. They found that for param- 
agnetic impurities the Hall conductivity vanishes. Our 
more general calculations confirm this result. However, 
we point to one important difference in its derivation. 



In both cases the dc-limit Kubo formula, where the con- 
ductivity is expressed via retarded and advanced Greens 
functions, has been employed to calculate the Hall con- 
ductivity. As was shown by Streda this version of 
the Kubo formula contains two parts: a^y a contribution 
from the Fermi surface and a^Jy a contribution from all 
states of the Fermi sea. The latter part is less known 
because it does not appear in the expression for the lon- 
gitudinal conductivity. Inouc et al}- calculated only a^y 
and indeed we find that for their choice of parameters the 
second part of the conductivity a^y vanishes, explaining 
the agreement with our results. In a more general analy- 
sis, beyond the limit of weak spin-orbit and Zeeman cou- 
plings, we find a non- vanishing a^y. Our work provides 
the missing estimate of cr^^ and extends the calculations 
of Inoue et ali^ 



Finally, the latest work on the subject is by Onoda et 
al}"^ The authors used the Keldysh technique, which they 
reformulated in a way appropriate for multiband prob- 
lems in a gauge invariant formalism. They also derived a 
self-consistent equation, which is the analog of the stan- 
dard quantum Boltzmann equation and solved it numer- 
ically. Unfortunately, lacking a full understanding of the 
details of the numerical procedure and the starting equa- 
tions being very formal within a non-chiral basis, a de- 
tailed discussion of their approach cannot be performed 
here. However, being devoted to the same model, the 
final results can be compared directly with the possible 
discrepancies arising from the different limits considered 
in the disorder distributions in which rii and the dis- 
order strength are two independent parameters in their 
calculations. Onoda et ali^ find a strong skew scatter- 
ing contribution of the order of esoyimpO'xx/W^ , where 
W is the inverse density of states. The skew scattering 
term changes sign at the point where the minority band 
becomes depleted, which they call the resonance point. 
The authors find also that the side-jump contribution is 
small in comparison with the intrinsic one. Our results 
confirm neither of those predictions. We find that for the 
Rashba model with randomly placed delta-function im- 
purities the leading part of the skew-scattering vanishes 
identically when the Fermi level is above this resonance 
point. Although skew scattering could still appear in 
higher order terms of the Born series, we expect these 
contributions to be small because they are of higher or- 
der in Vimp- On the other hand, Onoda et aS^ consider 
the limit of dilute impurities — > independently of the 
disorder strength Vimp which might be the origin for the 
discrepancies. Using the Keldysh formalism in the dis- 
order free basis we have been able to verify analytically 
our results. Further numerical analysis?''' of different lim- 
its will be necessary to settle the discrepancies with the 
results by Onoda et alJ^ 
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III. ANOMALOUS HALL CONDUCTIVITY OF 
THE 2DEG 

A. Model Hamiltonian 

We consider a spin-polarized two dimensional electron 
gas with Rashba spin-orbit interaction 



H = ^o-Q + oi{a^ky - aykrc) - ha^ -f V{v)<7f) (4) 

where m is the the effective in-plane mass of the quasi- 
particles, a the spin-orbit coupling parameter, h the ex- 
change field, and ai the 2x2 Pauli matrices. The eigenen- 
ergies of the clean system are 



Ek±^—± \k with Afc = V/i2+a2fc2 (5) 
Zm 

and are shown in Fig. [TJ The retarded Greens function 
of the clean system is: 



QiO)R ^ 



iO+] (7q + akyUx - akxO-y 



ha. 



V, 2m / 

G(")«ao + Gi°)^cT. + G(°)«a, -I- G(°)«a. , (6) 



with 



1 h 



Gr'' = ^(G?VGL°)) gW^ = -^-^(G 
2 2 Afc 



G 



(0)n 



2 A 



1 a/jj 

2 A/c 



and 



^:0^ 



(7) 

(Gf- gL"^) 
(8) 



The disorder potential V{r) in Eq. ^ is assumed as spin- 
independent. We consider the model of randomly located 
^-function scatterers: V(r) = ^ ■ Vi6{r — Hi) with ran- 
dom and strength distributions satisfying (Vi)dis — 0, 
{K^)d^s ^0 and (V^/)d.. = V-' ^ 0. This model is 

different from the standard white noise disorder model 
in which only the second order cumulant is nonzero; 
(l^k'kP)'iis ~ niVg where Ui is the impurity concentra- 
tion and other correlators are either zero or related to this 
correlator by Wick's theorem. The deviation from white 
noise in our model is quantified by Vi ^ and is neces- 
sary to capture part of the skew scattering contribution 
to the anomalous Hall effect. 

We calculate the self-energy using the Born approxi- 
mation: 



= -~i{T(To + VzCTz) 



(9) 




(b) 



FIG. 1: Single particle dispersion for small spin-orbit in- 
teraction akp/h — 0.2 (a) and large spin-orbit interaction 
akp/h = 2.0 (b). 



where i'± is related to the density of states at the Fermi 
levels of the two subbands 



dEk± 



dk 



(10) 



with 



A± = JK^ + a^kl = J\l + {a'^mf T a^m , (H) 



where Xp ~ \/h? + 2a^mep and 



k± ^ J2m [ er + a'^rriT \ >^F + (a'^m)'^ ] (12) 



are the Fermi momenta of the two subbands. 

Including the self-energy, the impurity averaged 
Greens function becomes: 



uj - + iTj (7o + akyUx - ak^Oy - [h + iTz)^^ 
Gq CTq + G^Gx + Gy Gy + Gf az ■ (13) 



By comparing this expression with Eq. ([6]) one observes 
that the impurity averaged Greens function can be ob- 
tained from the Greens function of the clean system by 
the following replacements: 

uj^uj + iT, h^h + iT^. (14) 

In the limit of small F^ one can therefore expand 



Afe ^ Vih + »r.)2 + a2fc2 ^Xk(l + i^] . (15) 



Using this approximation the impurity averaged Greens 
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function can also be written as: 
1 

2* 



= sin(?!)G'" = ^ 



with 



and 



(16) 



R 



1 aA;„A 



G? = 



cos0G^ = --^^2 



[Gl - G^) 



[Gf - G'^) 



2 A| + zr,/. (^+-^-^ 



1 



(17) 



(18) 






FIG. 2: Diagrammatic representation of the bare bubble (a), 
of the ladder vertex corrections (b) and of the skew scattering 
contribution (c). 



B. General expression for the anomalous Hall 
conductivity 

According to the Kubo-Streda formalism^° the ofF- 
diagonal conductivity can be written as: 



Jib) 



71 



where 

lib) 



2ttV 



a 



yx 



'yx 



AttV 



AttV 



TvivyG'^iepHG^ieF)) 

' Tr(«,G^(ef)i;,G«(eF) 
+VyG^{eF)vxG^{eF)) 

def{e)TT{vyG"{e)v, 



(19) 



(20) 



dG^ie) 



-Vy^^^VxG%)-VyG^ie)u 
dG^{e) 



de 
dG^je) 
de 



de 



VxG^ie)) 



Here, results from the electrons at the Fermi surface 
whereas a^^ denotes the contribution of all states of the 
Fermi sea. For cr^^^) and a^^ it is sufficient to calcu- 
late the bare bubble contribution in the weak scattering 
limit^" because vertex corrections are of higher order in 
the scattering rate F. Plugging in the Greens function of 
Eq. pB]) and using the velocity vertices 



The bare contribution of a" in the clean limit, i.e., for 
r+ = r_ = 0+ can be calculated by integration (see 
App. El and yields 



// = i_ 1 



where dG^^^ /de 



2a^mep + (a'^m)^ 



_^qR/A,, 



e(/i-eF)(23) 



has been used. Including 
the real scattering rates F+ and F_ does not lead to 
qualitatively different results but mainly causes a slight 
smearing. Thus we consider it as sufficient to focus on 
the clean limit contribution of <t^^ . 

For (T^(°) vertex corrections can be of similar magni- 
tude as the bare bubble and thus have to be considered 
carefully. In the weak scattering limit contributions of 
higher order impurity scattering vertices are small leaving 
only ladder type vertex corrections and the /{riiVQ) 
skew scattering contribution as the important termsJ^ 
Thus we decompose the following way: 



Ji^) 



_ ^ /(a),6 7(a)J 
^ yx yx ' ^ yx 



rHa) 



(24) 



I{a),b 



Vx 



m 



-o-Q - a<Ty 



one finds that cr^^''^ vanishes 

5^ 1 



a: 



Hb) _ 



yx 



4:TtV (27r)2 
-ia^G^G^ 



where ayx is the bare bubble contribution (Fig.l^^a)), 
ao + aa^ (21) CTyi"^'' the ladder vertex corrections (Fig. [IJb)), and 
CTyi"^''* the skew scattering contribution (Fig.^c)). With 
respect to the skew scattering contribution we have 
showni^ that only the diagrams with a single third order 
vertex (see Fig. ^c)) contribute to order Vi/{niVQ). In 
this diagram both vertices have to be renormalized by 
ladder vertex corrections. 



ia^G^G^) = { 



ia^G'^G^ 





(22) 
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1. Bare bubble 



and thus 



The calculation of the bare bubble contribution pro- 
ceeds as follows: 

= £// 9^TrKG-(..)..G-(e.)] 
dkk 



(1 - n^Vo^h)^ - [n^V^hY 
1 - TiiV^Ii iuiV^h 



(30) 



The ladder diagrams are therefore given by 



= 2m(2/3 - a/2) (25) 

where (for explicit evaluation of integrals Ii, I2, I3 and 
I4 see App. [B|) 



(26) 



/2 = 



1--)- 



i- 1 dfcfc(G?G^^-GfG^? 



27r 
1 

' 2^ 
1 



— / dk— g :g 



2m 



2m 



A+ 



dk 



r+A, 



r_A- 



A_ 



^. Ladder diagrams 

For the ladder terms a^i"^'' we sum the vertex correc- 
tions in front of the vertex as indicated in Fig. (Hb). 
Starting from the momentum integrated bare velocity 
vertex 



f f dkkd(f> ft, . ^A, , 
j J (27r)2 y.^F)'"^'-^ \^f) = 1^^^ + ^y^v ' 

with 

Ix = iih - ah) , jy ^ h - ah 
one finds for the renormalized vertex 

= Ixf^x + lyCTy 



(27) 



(28) 



(29) 



^^^G"{€F){jxTx+7yCTy)G^ieF) 
"JxCTx + lyf^y 

+niVQ {{h^x + ih^y)'yx + [h'^y - ih^x)cry) 



= -^^hy^x+lxTy) 



(31) 



n^Vp^ (27.7, (l-n^Fp^ Ji )+in,Vo^h (7^-7^)) 
{l-n,Vihr-{n^V,^h)^ 
In the weak scattering limit this reduces to 



riiVp {2'^xly(.l-niV^h)+iniV§hl^ 



(1 - n^Vihf 



Skew scattering 



(32) 



For skew scattering we consider only diagrams with 
a single third order impurity vertex and both external 
current vertices renormalized by ladder vertex corrections 
as indicated in Fig.^^c). In analogy to the renormalized 
Wa;-vertex in Eq. (|29p also the renormalized ii^-vertex can 
be calculated and expressed via F^^ and Fj, as 

(33) 



^ y^x 



^ x^y • 



Using these expressions the skew scattering diagram of 
Fig. IHc) yields 



27r 2 
2^~V^ 



^y"dfcfcTr[F,„G«(eF)F,^-f F„„F,^G^(eF)] 
Tr[-F^,^(F(To + ^zO-z)Tv^ 

+rt, F^^ (F(To + TzCTz)] 



e2 



■^^iTzT^'i'[{Ty(Jx + ^x0'y){0-z{rx0-X + TyCTy) 

Tyay)az) 



-(FxCr^r 



^ili 4F,(F2 



F" 



(34) 



From this expression it is evident that the skew scattering 
contribution vanishes as soon as F^ = implying that the 
lifetimes in both bands become equal since F_ — F+ = 
Tz{h/\^ -I- h/X+) vanishes for F^ = 0. Plugging in F^; 
and Fj, from Eq. ((30)) one findsi^ in the weak scattering 
limit, i.e., neglecting higher order impurity terms: 

s2 4lfF,72 



I(a),s ^ 
'yx 



2T:Vg {I -n,V;}hY 
2^n,Vo4jy_(3/i2 + A^)2 



(35) 
(36) 
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> 





with 



and thus 



T{a),b , 



T(a),l 



0, 



(42) 



i.e., the contribution of the bare bubble and the ladder 
diagrams cancel mutually. 



2. Only majority band occupied 



FIG. 3: Full vertex including ladder and skew scattering dia- 
grams. 



It can be shown easily that considering the weak scatter- 
ing limit of the full vertex shown in Fig. [3] yields exactly 
the same result as Eq. i.e., to order Vi/{niVQ) it 

reduces to the elementary skew scattering diagram de- 
picted in Fig. [DJc). 



C. Simple limits 

1. Both subbands occupied 

In the situation that both subbands are partially oc- 
cupied, i.e., ep > h, all contributions to the anomalous 
Hall conductivity vanish. For cr^^ this is immediately evi- 
dent from Eq. (|23|) . For the skew scattering contribution, 
which is proportional to Vz (see Eq. (|36|)). one observes 

easily that cjy^x^'^ = because (see Eq. ([9])) due 

to v+/\+ - iy_/A_ = (see Eq. CT). 

With respect to the bare bubble and ladder diagrams 
we will show in the following that they cancel mutually. 
For ep > h the integrals in Eq. (pS)) simplify to 



2 9 

a m ep 

2A|,r 



ihm 
2^ 



0, 



am , 
(37) 



and the bare momentum integrated vertices in Eq. 
are: 

a^mep\^ 



= a 1 



This gives for the bare bubble in Eq. ((25l 



I{a),b 



2tt 



^p 



(38) 



(39) 



For the ladder diagrams we need also 



-^n,V,^I, = (40) 
a 



yielding 



T{a).l _ 



2^ Xl 



Ix + 7y 



(41) 



In the opposite situation, where only the majority band 
is partially occupied, we have i/+ = and therefore F^ ^ 
0. In this case all terms contribute to the anomalous Hall 
conductivity. In the following we restrict our analysis to 
Fermi energies ep > —h, i.e., we disregard the region 
of very small Fermi energies, where the valley structure 
of the majority band becomes important (see Fig. [Hb)) 
and discuss the results in two simple limits: (i) small spin 
orbit interaction: akp <C h and (ii) small magnetization 
h <ti akp. 

In the limit of small spin-orbit interaction akp <^ h 
the sum of bare bubble and ladder vertex corrections be- 
comes 



e^iakpf 
27r IQhtp 



1 



(43) 



the contribution from the states of the full Fermi sea 



2 (akpf 
and the skew scattering term 



'yx 



27r ^epriiVa 



(44) 



(45) 



In the opposite limit of small exchange field h ^ akp, 
considering first a spin-orbit interaction still smaller than 
the Fermi energy akp <^ ep, we find for the sum of bare 
bubble and ladder vertex corrections 



a. 



I{a),b ,J{a),l 



yx 



3hep 



yx 



2tt {akpf 



(46) 



and for the contribution from the states of the full Fermi 
sea 

(47) 



4:71 \ akf 



and for the skew scattering term 



I(a)^s _ 
'yx 



2h£p 



2tt Vq n.iVoakf 



(48) 



In the same limit where the exchange field is small 
h <^ akp, but the spin-orbit interaction is now larger 
than the Fermi energy akp ^ ep we find for the sum of 
bare bubble and ladder vertex corrections 



a. 



I{a),b _^^I{a),l 



2he% 



2tt {akpY 



(49) 
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and for the contribution from the states of the full Fermi 



sea 



Sec. IIII C 11 Furthermore we observe that not only cr^^ 
but also the bare bubble and ladder vertex corrections 



/(a),b 



I(a).l 



TI 



e 

An 



1 



2he, 



and the for skew scattering term 



2tt V§ n,Vo 



(50) 



(51) 



D. Discussion 

We now discuss the full evaluation of the anomalous 
Hall conductivity in the limit of small spin orbit inter- 
action akp <C h and in the opposite limit of strong spin 
orbit interaction akp ^ h^ep- For the following discus- 
sion we will express all quantities in terms of the exchange 
field /i, which we define as /i = 1. Furthermore we will 
set m = 1, we choose Vi — Vq and use an impurity con- 
centration of Uj =0.1. 



cfyx " + cfy'x " (see Eq. P5|) ) are independent of impu- 
rity scattering. Both contributions are small: cr^^ con- 
tains a small prefactor of {akp/Kf' (see Eq. (01])) and 

small prefactor of {akp)"^ / [hep] (see 
Eq. (051) ). The skew scattering contribution, on the other 
hand, has a prefactor of akp/{niVo) which diverges for 
Vb ^ 0, i.e, l/r — > (see Eq. (05])) and therefore over- 
compensates the small prefactor oi akp/ep (see Eq. (05]) ) 
when the impurity potentials Vq becomes small enough. 
Thus for the parameters chosen in Fig. 0] the skew scat- 
tering term outweighs the other contributions by orders 
of magnitude and therefore the total anomalous Hall con- 
ductivity is almost identical to the skew scattering term. 
It increases quadratically with ep/h (see Eq. (0S])) and 
then vanishes suddenly for ep > h. 




Ef /h 



£f /h 



FIG. 4: Anomalous Hall conductivity for akp/h = 0.2 and 
an impurity concentration of rii = 0.1 plotted as a func- 
tion of ep/h (from right to left) and as a function of 1/r = 
rLimVo in units of h (from back to front) where upper left 
panel: total anomalous Hall conductivity (Eq. (|19|l ). up- 
per right panel: skew scattering contribution (Eq. (|36|l ). 
lower left panel: bare bubble plus ladder vertex corrections 
(Eq. (gg+Eq. lower right panel: cr" (Eq. AU 

conductivities are plotted in units of e^. 

In Fig. 0] we show the anomalous Hall conductivity for 
a small spin orbit interaction of akp/h — 0.2 as a func- 
tion of the Fermi energy ep/h and the scattering rate 
1/t = riiVQm for an impurity concentration of rii — 0.1. 
The upper left panel shows the total anomalous Hall con- 
ductivity, i.e., the sum of skew scattering (upper right 
panel), of bare bubble and ladder diagrams (lower left 
panel) and of the contribution from the whole Fermi sea 
(lower right panel). Obviously all contributions to the 
total conductivity vanish for ep > h, i.e., when both 
subbands are occupied which agrees with our analysis in 



£f /h 



-""IPIG. 5: Anomalous Hall conductivity for akp/h = 10.0 
and an impurity concentration of — 0.1 plotted as a 
function of ep/h (from right to left) and as a function of 
1/r = riimVo in units of h (from back to front) where upper 
left panel: total anomalous Hall conductivity (Eq. (US}), 
upper right panel: skew scattering contribution (Eq. (|36|l ). 
lower left panel: bare bubble plus ladder vertex corrections 
(Eq. ^+Eq. lower right panel: a" (Eq. AU 

conductivities are plotted in units of e^. 

Fig. [5] displays the anomalous Hall conductivity in a 
similar way as Fig.0]only for a large spin orbit interaction 

of akp/h = 10. A gain, (Tyx ^' -\- ayx ^' turns out to be 
independent of the impurity parameters and even smaller 
in magnitude as before because now it is suppressed by 
a small prefactor of {hep)/{akp)'^ (see Eq. (jl^ ). Analo- 
gously to the limit of small spin orbit interaction, the to- 
tal anomalous Hall conductivity is dominated by the skew 
scattering contribution, which contains no small prefac- 
tor and due to the factor of h/{niVo) grows rapidly for 
small impurity potentials Vo 0, i.e., 1/t ^ (see 
Eq. (|5T]) ). In the limit of large spin-orbit interaction 
akp ^ ep the skew scattering and thus the total anoma- 
lous Hall conductivity is independent of the Fermi energy 
ep ior ep < h (see Eq. (|5ip ) and then abruptly drops to 
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zero for ep > h. 



APPENDIX A: INTEGRATION OF a' 



IV. CONCLUSIONS 



In summary, we have investigated the anomalous Hall 
conductivity in a spin-polarized two-dimensional electron 
gas with Rashba spin-orbit interaction in the presence 
of pointlike potential impurities. Our calculations have 
been performed within diagrammatic perturbation the- 
ory based on the Kubo-Streda formula, an approach, 
which has previously been shown to yield equivalent re- 
sults to the semiclassical Boltzmann treatment j^i^^ 

Comparing our results with previous calculations we 
have been able to sort out contradictions existing in the 
literature. We have found that within the model Hamil- 
tonian considered all contributions to the anomalous Hall 
conductivity vanish as soon as the minority band be- 
comes partially filled, i.e., as soon as the Fermi energy be- 
comes larger than the internal Zeeman field. For smaller 
Fermi energies all contributions are finite with tr^^, the 
contribution from all states of the Fermi sea, being the 
smallest term at least in the limits of weak and of strong 
spin orbit interaction. The vertex corrections, which play 
the role of a side jump contribution, can be of similar 
magnitude as the intrinsic contribution and turn out to 
be independent of the impurity concentration and im- 
purity potential at least in the limits of small and of 
strong spin orbit interaction. In the weak scattering limit 
the dominant contribution results from skew scattering 
because due to its l/(niVb)-dependence it outweighs all 
other terms. Moreover, the intrinsic and the side jump 
terms contain higher orders of small prefactors than the 
skew scattering contribution. 
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Starting from the expression of a" in Eq. (I20p one 
obtains after angular integration: 



4^2 



- fdkk /°de/(e)^4Im [G^G«(G^ - 



(Al) 

Now performing the remaining integrals in the clean 
limit, i.e., using r+ = r_ = S, yields 



1 



47r 27r 
Im 



dkk defie)A—{Ek--Ek, 

J J-oo 



{e-Ek++i6Y{e-Ek-+i5Y 

Am 2ti J Afe Ek- - Eu+ 

{T:6{Ek+ - ep) + Tr5{Ek- - ep) 
2 

[ln(£;fc+ - ep- iS) 

~\n{Ek-~ep-i6)]}. 
Substituting E^^ — E^- — 2Afc and using 

/k 
dk^TT5{Ek± - ep) 



(A2) 



'^^k± \ u I 2 

TTTO 



■■7T6{Ek± - ep) 



■e(ep - E-. 



and 



A±|A± ± a^m 



dk—j hi{Ek± — ep — iS) = 



(A3) 



\n{Ek± - ep - iS) 



dk 



1 



dEk± 



a^Xk Ek± — ep — iS dk 



J 

(A4) 



and 



ln{Ek+ - ep - iS) 



a'Xk 
-^eih-ep) 



lii{Ek- - ep - i6) 



a simplifies to 



= - — h 



1 1 



1 1 



47r V ~ orm m A+ + a^m 



--e{h-ep) 
h 



J 

(A5) 



Q{ep-h) 



47r \ i//i2 + 2a'^mep + (a^m)^ 
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APPENDIX B: INTEGRALS IN THE WEAK 
SCATTERING LIMIT 

In the weak scattering limit (F, small) the integrals 
over two Greens functions simplify to: 

^ I dkkmGl{k)Gl{k) (Bl) 

If 1 1 

= — / dkkfik) — 

27r7 •'^ ' eF-Ek++iV+eF-Ek+-iV+ 

2r+ 

^ / dkkf{k)G^{k)G^{k) « ^'l^''' analogously 

ZTT J 2r_ 



and 



i- j dkkf{k)Gl{k)G^{k) (B2) 

If 1 1 
= — / dkkf{k) — 



ilT5{€F - Ek-) 

1 



cf - Ek 

If 1 
- / dkkfik) 



efe — Afe EF — Efe + Afe 
1 



+ - / dkkfik) [5{tF-eu + Xk)- . 
2 ./ \ - Efe - Afe 



-(5(eF - Efe - Afe) 



CF — Efe + Afe 



yielding 



i- 1 dkkf{k){Gl{k)G^{k)-G^{k)Gl{k)) 
u-f{k{Ek-))5{eF-Ek-) 



i I dEk- 



i / dEk+ 



£f - Ek- - 2Afe(Ffc_) 
u+fik{Ek+))5ieF - Ek+ 
e-F — Ek+ + 2Xk{Ek+) 



i ( v+f{k+) ^ v-f{k-) 



2V A, ■ A_ 

Now we find for the integrals /i, I2, I3 and I4 in the 



weak scattering limit: 

^1 = ^/ dkk{G^G^-G^G^) (B4) 

= \^Jdkk {G^G^ + G^G^ + G^G^ + G^G^ 
Xlih' + Tl) 



{GlGl + G'^G^ - GlG^ - G'^Gl)^ 

/2 = ^ y dkk {G^Gf - G^Gi) (B5) 

+ iT,{h^-Xl){G^Gi-G^G^)) 
Jdkk^^{hXU-G^G^ + G^Gi) 

+iT4h^ - XDiG'lG^ - G^G^)) 
I (u+h v-h y+{h'^-X\) v-jh'^-Xl.) 

4VA+ ^ +r+ x\ r_ A3 



=---[ 

2 27ry 



n/F X2 



'2m A^+r2/i2 



(G^G^ + G!!G^ (B6) 
-G^G^ - G^G^) 



.1^ / dkk^^iG^Gi + G^G^) 
22tt J 2m A2 ^ + + " "-^ 



(B7) 



2 27ry 



1 1 



dkk 



fc^ aAfe 
2^A^ + r2/i2 



(A|(G^G^ - G^G^) 



+zr,/i(G^G^ - G^G^)) 



.1^ / dfcfc^ii(G^G^-G«G^) 
2 27ry 2mAfe^ + + ~ ~' 



■^"'^^'TvA;-Ta_ 
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